arXiv:1501.02748v4 [cond-mat.mes-hall] 22 Oct 2015 


Majoranas with and without a ‘character’: hybridization, braiding 

and Major ana number 

N. Sedlmayr/’ * M. Guigou,^ P. Simon,^ and C. Bena^’^ 

^Institut de Physique Theorique, CEA/Saclay, Orme des Merisiers, 91190 Gif-sur-Yvette Cedex, France 
^Laboratoire de Physique des Solides, UMR 8502, Bat. 510, 91405 Orsay Cedex, France 

(Dated: October 23, 2015) 

In this paper we demonstrate under what conditions a pseudo-spin degree of freedom or charac¬ 
ter can be ascribed to Majorana bound states (MBS). These exotic states can be created at the 
boundaries of non-interacting systems, corresponding to D, Dill and BDI in the usual classification 
scheme, and we focus on one dimension. We have found that such a character is directly related to 
the class of the topological superconductor and its description by a Z, rather than a Z 2 , invariant 
which corresponds to the BDI class. We have also found that the Dill case with mirror symmetry, 
which supports multiple MBS, is in fact equivalent to the BDI class with an additional time-reversal 
symmetry. In all cases where a character can be given to the Majorana states we show how to 
construct the appropriate local operator explicitly with various examples. We also examine the 
consequences of the Majorana character by considering possible hybridization of MBS brought into 
proximity and find that two MBS with the same character do not hybridize. Finally, we show that 
having this character or not has no consequence on the braiding properties of MBS. 

PACS numbers: 73.20.-r, 73.63.Nm, 74.78.Fk 


I. INTRODUCTION 

In his seminal 2001 paper, Kitaev introduced a simple 
toy model of a Id topological p-wave superconductor of 
length L.^ Local Majorana operators were introduced by 
decomposing each local annihilation and creation elec¬ 
tronic operators Cj and c| at a given site j into a pair of 
local Majorana operators defined by ^a,j = (cj -|- c])/2, 

Ib.j = {cj — cj)/2i. The topological phase is charac¬ 
terized by two localized Majorana fermions at both ex¬ 
tremities of the chain. When the pairing and hopping 
amplitudes are equal, the Majorana operators 70,1 and 
"fb.L decouple from the Hamiltonian. These are the zero- 
energy modes characterized by Majorana operators ja.i 
and 7 b,L. Although in this decomposition, the indices a, b 
may appear somehow arbitrary (note also that a and b 
can be exchanged), they are related to the trivial feature 
that a Majorana operator is always built from two differ¬ 
ent fermionic operators in any realistic condensed mat¬ 
ter Hamiltonian. One may wonder whether - and under 
which conditions - we can associate to this label a kind 
of pseudo-spin degree of freedom or ‘character’. In the 
Kitaev chain, if a Majorana fermion of type ‘a’ is local¬ 
ized at one end of the chain, then necessarily a Majorana 
fermion of type ‘b’ is localized at the other extremity. 
In this paper, we will be focusing on non-interacting Id 
systems. 

In order to attribute some physical meaning to this 
‘character’, we are looking for a local operator V which 

(i) anti-commutes with the Hamiltonian H, {H,V} = 0, 

(ii) obeys = 1 , and (iii) is unitary and thus distinct 
from the particle-hole symmetry operator. Such an op¬ 
erator is an example of a chiral symmetry.^ The first 
two properties ensure that the Majorana bound states 
which can appear can be characterized as eigenstates of 


V with two possible eigenvalues defining the aforemen¬ 
tioned ‘character’ of the Majorana fermion. Moreover, 
it follows that if the Hamiltonian has reflection symme¬ 
try, then these Majorana bound states will be localized 
at the two boundaries. Any continuous change to the 
Hamiltonian which does not close the gap and leaves the 
symmetry intact will preserve this. Thus we always have 
localized Majoranas with a definite character, provided 
the above conditions are satisfied. 

Since we focus on Id particle-hole symmetric non¬ 
interacting systems, they belong to one of the following 
three classes: BDI, D or Dili depending on the presence 
or absence of different anti-unitary symmetries.^ These 
symmetries are referred to as time-reversal symmetries 
and at the single particle level satisfy T± = ±1 for 
BDI/DHI and are absent for D. In systems belonging 
to the topological class BDI, the symmetry associated to 
the operator V is the composite of the time-reversal and 
particle-hole symmetries and is also known as the chiral 
symmetry, the existence of which is already well-known.^ 
We prove the equivalence between the chiral symmetry 
and our character operator for this class. The character 
however can also be viewed as a local property of the 
MBS, a property we will extensively exploit to analyze 
various examples. Each ‘type’ of Majorana is localized 
on a different boundary of the system. In the class D, one 
no longer has the symmetry which allows this classifica¬ 
tion. In Id the existence of this character can be directly 
related to the change in topological invariant from Z for 
BDI to Z 2 for class D as the possibility of having many 
Majorana states at the end of a wire^ requires them to 
have a well-defined character. This is because, as we will 
demonstrate, Majorana states without such a character 
generically destroy each other on contact, recombining to 
form finite energy Dirac particles. 

Dili topological superconductors are time-reversal in- 
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variant and have a Z 2 invariant.Thus generically they 
behave like class D but with Kramer’s pairs of Majorana 
states. However in addition to time-reversal symmetry 
they can also still have a so called mirror symmetry j\4 
where = —1.® We show that when [H,M] = 0 one 
also has {H,V} = 0 and a character can be defined. In 
these cases the same conclusions for the BDI class can be 
applied to the Kramer’s pairs of Majorana bound states 
and one can have many Majorana states. When this mir¬ 
ror symmetry is broken, then {i7, V} 7 ^ 0 and the Dill 
class behaves like class D, which are both described by a 
Z 2 invariant. The possible existence of the many Majo¬ 
rana states in the presence of the mirror symmetry indi¬ 
cates that it should have a Z invariant,® and we show that 
this case is in fact equivalent to the BDI class with an ad¬ 
ditional time-reversal symmetry which does not change 
the symmetry class. Note that we are not arguing simply 
that the low energy sector is BDI, as can arise in some D 
systems,^ but rather that in this case the full symmetry 
of the Hamiltonian is BDI, and not Dili, without further 
approximation. It has already been noted that multiple 
particle-hole symmetries or time reversal symmetries of 
the same type do not alter the classification scheme as the 
multiple of, for example, two time-reversal operators is a 
unitary operator.^ However this does not in itself decide 
which class a system with two different types of time re¬ 
versal symmetries belongs to, and consequently what its 
topological invariant is. The simultaneous presence of 
two different time-reversal symmetries is still ambiguous 
in the Schnyder et al. classification scheme,^’® and we ar¬ 
gue here that this class should be properly thought of as 
a subset of BDI. We prove that although the Dili class 
also possesses a chiral symmetry this can not be related 
to a character. 

The interest in Majorana bound states is par¬ 
tially fuelled by their possible application to quantum 
computing.^’® To this end the non-abelian braiding prop¬ 
erties of the Majoranas are key.®°“^^ We consider the im¬ 
plications for the braiding properties for Majorana states 
which have no well defined character, where the usual 
mapping to a Kitaev chain is modified. Majorana states 
with definite character can be shown to either hybridize 
to finite energy states, or to be protected and remain at 
zero energy. Which of these occurs depends on whether 
the Majoranas have the same character, when they are 
therefore protected, or a different character when they 
will hybridize. Majorana states which exist without char¬ 
acter are much less robust and their mutual interaction 
will, unless one can carefully tune the appropriate over¬ 
laps to zero, destroy the zero energy Majorana states. 
The resulting finite energy states remain localized at the 
edges, but do not retain the Majorana anti-commutation 
relations which occur uniquely at zero energy. We go on 
to demonstrate that the character or absence of it has no 
known effects for the braiding properties of the Majorana 
states. Indeed, as one wishes for the states which are to 
be braided to hybridize between themselves, thus avoid¬ 
ing problems caused by degeneracies at zero energy, the 


case where the full system allows a well defined character 
needs to be avoided. 

The paper is organized as follows. In Sec. H, we exam¬ 
ine under which symmetry conditions we can build the 
Majorana character. In Sec. HI we demonstrate the con¬ 
sequences of the Majorana character by considering pos¬ 
sible hybridization and destruction of Majorana bound 
states brought into proximity. In Sec. IV we consider 
the implications of Majorana character or its absence on 
their braiding properties. We then summarize our results 
and offer some perspectives in Sec. V. 

II. CHARACTERIZING MAJORANA STATES 

Before characterizing the Majorana states, let us in¬ 
troduce our notations and conventions. 


A. Particle-hole symmetry operator and notations 

We will use the Nambu basis throughout the paper: 
'k] = -Cj-f-}, where annihilates (cre¬ 

ates) a particle of spin cr at a site j. The corresponding 
wavefunction is ipj = {uj^,Uji,Vji,Vj^}. We will also 
use T to denote the Pauli matrices in the particle-hole 
subspace and a as the Pauli matrices in the spin sub¬ 
space. We confine ourselves to discussing only generic 
Bogoliubov-de-Gennes Hamiltonians, H, which are anti¬ 
symmetric under the particle-hole transformation C = 
e^^a-yryK, i.e. {H,C} = 0 and = 1. Here K is 
the complex-conjugation operator, and C is an arbitrary 
phase which, without loss of generality, we will set to 
zero for convenience. Therefore the Hamiltonians have 
spectra (±ei, ± 62 ,..., )j where d is the Hilbert space 
dimension. 

Now consider two eigenstates of the Hamiltonian ^±1 = 
^(±ei) such that {H =F ei)^±i = 0. Let 

) ( 1 ) 

JO- 

then by applying the particle-hole transformation 

C-i = ■ (2) 

JO- 

Here a =t,4. should be taken as cr = ±1 whenever a 
numerical value is required. If the system is in a topo¬ 
logically non-trivial phase then Majorana states will be 
present at the boundaries.In the thermodynamic 
limit, system size L —>■ oo, the energy of these states 
goes to zero exponentially, limi_>oo Ci ~ where 

£ is a constant. It is then always possible to construct 
two Majorana states out of these degenerate solutions. 
These Majorana states can be localized at each bound¬ 
ary, though naturally other bases are possible. As a 
Majorana is its own anti-particle then a Majorana must 
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be an eigenstate of C with an eigenvalue of magnitude 
1, C I 7 ) = e*^| 7 ). Once again we can always tune the 
phase to be C = 0 by changing the irrelevant over¬ 
all phase of the eigenstates. Therefore if we calculate 
the expectation value, with 7 ^ = = 

then we find 

I < 7 IC I 7 > I = |2 ^ ^ crUjo-Ujcrl = |2 'y ' cr|ujo-| 1 = 1- 

i.o- 

( 3 ) 

This is equivalent to imposing the appropriate anti¬ 
commutation relations on the Majorana states. All non¬ 
zero energy eigenstates must satisfy 

<^n|C|e„ > = <^„|C-„ >=0. (4) 


B. General symmetry considerations 

The systems we consider generically possess the prop¬ 
erty that {i?, C} = 0 with = 1. We assume that we 
also have some V which satisfies {7J, V} = 0 for a given 
Hamiltonian H with = 1. Although the condition 
I V P = 1 would be sufficient, as V and C must commute 
the stronger condition = 1 is required. Additionally 
the operator V should commute with all other symmetry 
operators of the Hamiltonian. As the operators C and 
V are assumed to be distinct and we have C A = X* C 
where A is a complex number and A* the complex con¬ 
jugate, V must satisfy VA = AV. Now let us suppose 
we have two zero energy states ^± 1 , where as before for 
a finite system zero energy really means exponentially 
small in the system size. These states are related by 
CC±i = Cti and VC±i = however Cif±i = -i^^i 
and Vi^ii = iC=Fi- Constructing the orthogonal Majo¬ 
rana states 71 = (^1 -I- and 72 = i(^i - 

we immediately have that C 71.2 = 71,2 but that V 71 = 
7 i and V 72 = — 72 - The two eigenvalues ±1 define the 
Majorana character. 

There are three symmetry classes which we must con¬ 
sider. Firstly we have class BDI and class Dili which 
have a time-reversal invariance operator, T±, satisfying 
[H,T±] = 0 and T± = ±1 respectively. Then we have 
class D which has no time-reversal symmetry. 


1. Class BDI 

For the BDI class, the anti-commutation relation 
{H,V} = 0 and = 1 are ensured for V defined by 
V = CT+. V is sometimes called the chiral or sub-lattice 
symmetry. This symmetry allows one to assign a chiral 
character to the Majoranas as we shall see, and this in 
turn allows many Majorana states to exist at the end of 
a wire without recombining. This is in agreement with 
such a system having a Z invariant. In the simplest cases 
T + = K the complex conjugation operator and many 


other BDI cases can be transformed to this one by a 
global transformation, see App. A. 


2. Class D 

For class D, we find that it is not possible to assign a 
character to the Majorana states and only a single Ma¬ 
jorana can exists at one end of a topological wire. Such 
a system has a Z 2 invariant. 

To prove this, let us assume that we can define an oper¬ 
ator V that satisfies our requirements, namely: {H,V} = 
[C,V] = 0 and = 1. Then we can define an opera¬ 
tor T = CV. We see that such operator T must obey 
[H,T] = 0 and = 1. Therefore this means we are 
in the BDI class and not in the D class as originally as¬ 
sumed. Therefore this proves that class D has no charac¬ 
ter. Physically this makes sense if we regard the Majo¬ 
rana character as allowing the coexistence of other Ma¬ 
jorana bound states. As this is impossible in class D in 
Id due to the Z 2 topological invariant, there cannot be 
any character in class D. 


3. Class Dill 

For class DIB we shall consider two cases. Firstly if 
we do not have mirror symmetry then the system is in 
the same situation as for the D class and we still have 
no character. The same proof as for the D case can be 

given. We note that an operator defined by V = C T- 
- 2 

can not be used as V = — 1. 

The second DIB case corresponds to the Hamiltonian 
having an extra-symmetry, a so-called mirror symmetry 
introduced in Ref. 5. Mirror symmetry is defined by a 
unitary operator A4 which commutes with the Hamilto¬ 
nian, C and T- and also has = —1.® This is in fact 
a Hamiltonian with the symmetry [H, T+] = 0 in addi¬ 
tion to the physical time-reversal invariance of the DIB 
class [H,T-\ = 0. That mirror symmetry and simulta¬ 
neous T± symmetries are equivalent one can easily see 
by the relation T+ = T- M .. It follows directly from the 
properties of T_ and M that [H,T+\ = 0 and T+ = 1, 
which are of course the definition of T+. Similarly, if one 
starts with T+ and T-, it then naturally follows that 
one can construct an operator M where \H, M] = 0 and 
= —1. The operator V remains the same as for the 
BDI class and all properties of these states transfer to 
the Kramer’s pairs. As before V —CT+. 

This demonstration shows that a DIB system with a 
mirror symmetry should be thought of as a member of 
the BDI class instead. As multiple Majorana states are 
possible, it must be characterized by a Z and not a Z 2 
invariant, as was recognised early in Ref. 5. 
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FIG. 1. (Color online) Schematics of a superconducting wire 
a) in the presence of a perpendicular homogeneous magnetic 
field B and spin-orbit coupling a and (3, b) in the presence of 
a rotating magnetic field B. 


III. CONSEQUENCES OE THE CHIRAL 
MAJORANA CHARACTER ON 
HYBRIDIZATION 

To investigate the robustness of the Majorana states 
at the ends of a wire we will consider a ladder made by 
two coupled wires, with v = 1,2, described by the 
Hamiltonian 


H = 


Hi H,\ 

Hi H2 ) ■ 


(5) 


He is the coupling between the wires. As a first step we 
diagonalize H 12 using 


N = 


Ni 0 \ 

0 N 2 ) 


( 6 ) 


where 2 -^i, 2 A^i ,2 are diagonal. For we introduce a 
general one dimensional tight-binding Hamiltonian on L 
sites 


A, = Ho + Hb + + H, + Hd , (7) 

There are already several possible experimental realiza¬ 
tions of Majorana fermions, though definitive proof is still 
lacking. The models introduced here are sufficient to 
capture the essential physics of these various systems. In 
Fig. 1 schematics of the wire systems are presented. 

Suppressing for the moment the wire index /i, the hop¬ 
ping and chemical potential are 

^0 = E Km - t)r^] + H.c.j . 

( 8 ) 

We consider two different forms of superconductivity in¬ 
duced by proximity to a substrate. Firstly we have s- 
wave superconductivity, for the BDI and D topological 
superconductors, 


= (9) 

3 

For Dill we consider d-wave superconducting pairing 
induced by proximity with a dx^-y^ superconducting 


substrate" 






H.c. 


( 10 ) 


In both cases we have defined the pairing as 

. Tk A H . U: 


* S.d ^ A S.d _ 

2 ^ 2 




V 


( 11 ) 


t is the nearest neighbor hopping, fi a chemical potential, 
and the magnitude of the superconducting pairing 
assumed to be induced by a proximity effect. The super¬ 
conducting phase is Kj. Additionally there is a magnetic 
field of strength B, 


= ( 12 ) 

3 


which can locally vary its orientation in a direction given 
by the unit vector Such a magnetic field could 

for example be created by considering local magnetic 
moments of magnetic adatoms on a superconducting 
surface.For class Dili we must have 5 = 0 for time- 
reversal invariance. 

Finally we can have spin-orbit coupling 


Hso = -J2l + H.c.j . (13) 


Both Rashba and Dresselhaus terms are present, given 
in a low-energy approximation by acr^ and f3cr^ respec¬ 
tively. In this approximate form there is no difference 
between the effects of Rashba and Dresselhaus couplings 
and a global spin rotation around the z direction suffices 
to align the spin-orbit coupling along, for example, cr^. 
Note that in this model all terms except for the magnetic 
field and spin-orbit coupling are invariant under a global 
spin rotation. We set t = h = a = 1, where a is the 
lattice spacing. 

The inhomogeneous magnetic field, which can vary its 
direction arbitrarily, can be described by a unit vector 

fij = (cos Dj sin ipj, sin Dj sin ipj, cos (pj), (14) 

though we will generally focus on examples where 

Pj = 2'Kk,p{j — l)andi?_,- = 2TTk^{j — 1). (15) 

We are interested in two different cases with either spin- 
orbit coupling or inhomogeneous magnetism, though the 
composite can also be considered. 

If both wires are assumed to be in the topological 
regime then the matrix elements of n\HcN 2 which cou¬ 
ple the topological edge states to the bulk states are 
small, and the problem approximately decouples into a 
bulk and edge part. We are only interested in the Majo¬ 
rana states and so we focus on the effective Hamiltonian 


( -ei 0 M _M_+ \ 

0 ei M+_ M++ 
M*_ M*_ -62 0 

MI+ m;+ 0 62 / 


(16) 
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The matrix elements are given by 

Mab = {^{aei)\Hc\^{he2)) 


(17) 


- AM \M. 


where Hi 2 i{aei^ 2 ) = aei, 2 '?(«ei, 2 ) with a = ±1 and 
limL_>oo Cl ,2 —t 0. We will now explicitly use this last 
property, setting ei ^2 ~ 0. The eigenvalues of Eq. (16) 
are then 

f++|V|M+_|' + |M_+|V|M__|2) (18) 

+ \M++M__ = 0. 

Eq. (16) can also be written in the Majorana state ba¬ 
sis, which appears the same after a substitution of the 
eigenfunctions. We use both formulations in the follow¬ 
ing arguments. When all four matrix elements vanish, 
then all four Majorana states survive. When only the 
0(A°) term vanishes then only two Majorana states sur¬ 
vive, this is a finely tuned example and is not usually 
seen. The extension to cases where there are more than 
two zero energy states in the wires is straightforward and 
we do not give it explicitly. 

The overlap can be written for eigenstates where 


II 

-H 



> 

II 


jlcrcr' 

2 -^ 

2 


+ + c>) 




^ja' ^jo 


(19) 


We have used a coupling of the form 


( 20 ) 


which obeys particle-hole symmetry and has spin struc¬ 
ture t'. In order for this matrix element to vanish it is 
sufficient that there is the same chiral Majorana char¬ 
acter on any pair of two coupled sites belonging to the 
two wires and that t' does not break the chiral symme¬ 
try. This is one of the most robust ways to make the 
overlap vanish, though it is not the only possible way. 
This is clearest for the case where 5^=0 and Kj = 0, 
then the eigenstates of the Hamiltonian are real and the 
Majoranas are purely real at one end of the wire and 
imaginary at the other. Therefore all = 0 provided 
is short ranged and does not extend along the en¬ 
tire system. This argument extends directly to all planar 
spin textures and homogeneous superconducting phases, 
as the required transformations do not alter Mab- 


A. Coupling BDI wires 

We can now consider a variety of different scenar¬ 
ios starting with coupling BDI wires. First, by taking 
fij = {0, 0,1} and = 0, with ax = ±02 > 0 we can 


flip the chiral character of Majoranas in the wire 2 by 
changing this sign. Note that flipping a —> —a is not the 
only way of reversing the Majorana character, shifting 
from p, =—>■ 2 — /i has the same effect. By varying the 
coupling between the wires, t' we find the conditions of 
destroying the Majorana states. Second, by coupling the 
two wires at one single point, we can consider two wires 
of different types joined end to end. By calculating the 
hybridization of Majorana states in disconnected wires 
when they are brought together, it is possible to predict 
under what conditions Majoranas destroy each other or 
co-exist, and to relate this to the Majorana’s character. 
The same situation exists if we have a^, = /3,^ = 0 and 
a varying magnetic field confined to a plane, as can be 
seen in 2 d arrays. 

If {i7, V} 7 ^ 0 then only for fine tuned examples can we 
have Mab = 0- Generically the Majoranas will hybridize 
and gain a finite energy. This follows from the obser¬ 
vation that M^l is no longer fixed to zero and therefore 
Mab = 0 requires either a careful cancellation in the sum 
over position or a tuning of the coupling t' ^. In this case 
the full system is in the D class, even though each wire 
separately is BDI and the coupling can be trivial. 

Although it is true to say that Majorana states with 
the same character will not hybridize, if they have oppo¬ 
site chiral character the overlap emphcan still be zero, 
again this requires some fine tuning. For example if 
/ii = 0 and fj ,2 = ‘2 with all other quantities being equal 
between the wires, then the Majorana states of the inde¬ 
pendent wires will be orientated in opposite ways. Conse¬ 
quently in one wire the state V = 1 will be on the left and 
in the other wire it is the V = — I state that will be on the 
left. Provided = t'Sij we still have Mab = 0. A small 
amount of disorder is however enough to couple the states 
giving Mab 7 ^ 0. This is an example of weakly protected 
topological states in a topologically trivial phase, referred 
to elsewhere as hidden symmetry.In fact if {H, V} 7 ^ 0 
then Mab = 0 due to two possible reasons. Firstly there 
can be an emergent symmetry in the low energy sector in 
which case effectively we recover {H,V} = 0 for an ap¬ 
proximate low energy Hamiltonian H, and secondly we 
can have some form of fine tuning of the coupling or weak 
topology. 

The case where we have both [H,T-] = 0 and 
[H, T+] = 0 is similar to the simple BDI case where 
the states now exist in Kramer’s pairs protected by the 
symmetry. For example a system with d-wave pair¬ 
ing, Rashba spin-orbit coupling, and no magnetic field, 
i.e. Hi, with B = A® = 0 has both [iJ, T+] = 0 
and [B, T+] = 0. It can support multiple Kramer’s 
pairs of Majorana states at each end with a character 
V = CT+ = cT^Ty. Breaking T- takes us back to the 
usual BDI case with no Kramer’s pairs. Breaking T + de¬ 
stroys the symmetry protecting the multiple Majoranas 
and we can have only a single Kramer’s pair at each end 
of the wire. 

We can summarize the findings in the following way. 
For wires in the BDI class then the Majorana states are 
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protected provided the ladder has the same symmetry 
and the wires are weakly coupled. Strong coupling can 
drive the system through a topological phase transition 
back to a topologically trivial regime. This is not possible 
with weak coupling as the topological phase transition 
requires a closing of the bulk gap. If the symmetry is 
broken then all the Majorana states are normally gapped, 
though if the ladder remains topologically nontrivial then 
a single Majorana at each end survives. 


B. Coupling D and Dill wires 


If we start with wires with D symmetry then, except 
for very carefully tuned examples, the Majorana states 
are unprotected and become gapped. The Dill case is 
essentially similar to the D case. We remind the reader 
that the Dill symmetry with mirror symmetry is consid¬ 
ered as a subset of BDI (see Sec. II). 

If we have two wires with slowly varying supercon¬ 
ducting phases Kj, where the variation is confined to the 
bulk of the wires then we are in the D class. Now as we 
can locally define a character at a single end of the lad¬ 
der we should be able to have multiple Majorana states, 
although we are in D. In fact such a situation has an ef¬ 
fective low energy theory in the BDI class and the Majo- 
ranas can only survive to the extent that one can neglect 
the corrections to this effective model which lift it out of 
the BDI class. 


C. Coupling Kitaev chains 


If we consider coupling Kitaev chains rather than full 
spinfull Hamiltonians then we can consider other pos¬ 
sibilities. One can define a time-reversal operator which 
operates in the sublattice space of the ladder, T- = X^K 
where are Pauli matrices operating in the ladder 

subspace. It is then possible to construct Dill ladders 
from D wires and BDI ladders from BDI wires, addition¬ 
ally we can impose the mirror symmetry in this second 
case. With the mirror symmetry the many possible Ma¬ 
jorana states which can be formed at the end of a wire 
are protected. For class D the states are in fact protected 
by the time-reversal invariance of the Dill ladder, if this 
is broken then they will become gapped as predicted by 
the absence of a well defined character. 

We now give explicit examples. The full system Hamil¬ 
tonian can be written as 

The wires, where we have included longer range hopping 


and superconductivity , are 

3 


( 22 ) 


< 

0 e'^^o ^ 

- 

f'J \ 

— 0 J 





where = {(3^pCvj}, where annihilates (creates) a 
spinless particle at a site j on wire v. We also introduce 
This chain can support two Majorana 
bound states when its symmetry class is D, and many 
Majorana bound states when its symmetry class is BDI. 

If we consider two BDI wires with the symmetry 
= 0 where in this case V = r^, then provided 
{H, V} = 0 all Majorana states are protected and the full 
system is still BDI. For example 


= (23) 

j 


A coupling which breaks {H, V} = 0, or a combination 
for different BDI wires which breaks this symmetry, ends 
up in a D ladder with no Majorana states surviving. Note 
we are by construction in the topologically trivial regime 
of the D ladder. 

We can also construct a coupling which gives us the 
effective time-reversal invariance, [H,T-] = 0, for T- = 
X^K, given by 

H^ = A '. (24) 

j 

If the wires we couple are in the BDI class then this 
extra symmetry corresponds to mirror symmetry in the 
ladder Hamiltonian, and we have many Kramer’s pairs of 
Majorana states. If we couple D wires then we have no 
mirror symmetry and only a Kramer’s pair of Majorana 
states is possible, in agreement with the symmetry class. 
Note that in this case the Majorana bound states of the 
uncoupled D chains survive as the Kramer’s pair in the 
Dili ladder, they are protected by the time-reversal in¬ 
variance. Longer range coupling between the wires does 
not change these conclusions. 

We have extensively checked all the results of the anal¬ 
ysis in this section, which is based on the hybridization 
of the zero energy states inside the gap, by performing 
an exact scaling analysis on the low energy states as a 
function of system size. Numerically solving the tight 
binding Hamiltonian Eq. (20) we can unambiguously test 
whether the Majorana states survive. In all cases, pro¬ 
vided inter-wire coupling is weak t' ^ t/5, the results are 
all verified. 


IV. BRAIDING CHARACTERLESS 
MAJORANAS 

To understand if the characterless Majoranas 
{{H,V} 7 ^ 0 ) would braid differently we consider a 
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generic example of such states and map the correspond¬ 
ing model to the Kitaev model. In the limit of large 
magnetic fields it is possible to map the spinfull problem 
to a simpler low-energy model similar to a Kitaev 
chain. First we consider the continuum limit of our 
model for slowly varying magnetic fields. After a gauge 
transformation which fixes the local spin direction to 
be parallel to the magnetic field, see App. B, we have a 
homogeneous Hamiltonian with a vector potential. The 
vector potential is defined as A(r) = U\r)'S/U{r). This 
can be understood as the term which appears in the 
lowest order expansion of the effective spin-orbit term 
« <5 • A{Ri) where Ri is the real space position of the 
i site and once again Sij is the nearest-neighbour hopping 
vector from j to i. Naturally, for the one-dimensional 
systems considered here, such terms are non-zero only 
along the x-direction. We find 


to characterless Majoranas ({77, V} 7 ^ 0 ). The effective 
parameters become 


|Aeff| = 


sin(27rA:^x) 


-I- sin^( 27 rfc,^x)] , and 

(/left = tan“^ . (30) 


However such a phase variation plays no role in the braid¬ 
ing properties of the states, and can be deformed contin¬ 
uously without affecting the results, provided the gap is 
not closed.Thus we would expect that at a junction 
of three wires, the braiding operations with such char¬ 
acterless Majorana states remain the same as those pre¬ 
dicted for the simpler ({77, V} = 0) models. Moreover, 
our observation re-enforces the fact that the braiding of 
the Majorana states should not be affected by magnetic 
impurities. 


Ax {Rj) 


1 

2 


[dxipcr^ + Ox'S (cos (fcr^ — sin (pcr^)] ■ 


(25) 


V. CONCLUSIONS 


The Hamiltonian, Eq. (7) with a = P = = kj = 0 

and A®. = A, then becomes 

H f dxpjlx [e- f^ + Ba^^] ijj^^x 

a 

-A J dxtpi^xi^^^x + H.c. (26) 

-I- J dxV’t.x^ [dx^ - i^j^r^sint^] dx4>i,x + H.c.. 

iequiv — dxxli‘^'m) and the diagonal contribution to the 
gauge potential has been neglected. At low energy, and 
for sufficiently small spin-orbit terms this can be mapped 
to a spinless model via 

( 27 ) 

i’l,X ~ 'P'x 

resulting in the effective single band model 

B ^ J dxtpl [e - fj,es] i^x 

- J dx| Aeff|e““^“f‘'V'd;'9a:V’x + H.C. . (28) 

The effective parameters are 


Meff // T -S , 

|Aeff| = T (^2 ^ 


peS = tan 


dx'&shlip 

dxV 


(29) 


We have an effective Kitaev model with a spatially vary¬ 
ing superconducting pairing and phase. 

As an example we consider a combination of a process¬ 
ing and rotating field, given by Eq. (15), which gives rise 


We have investigated the possibility of characterizing 
the Majorana bound states which occur in topological 
superconductors. We find two distinct cases. On the one 
hand, which we call case (A), we have all systems with 
a time-reversal symmetry whose anti-unitary operator 
squares to 1. This is the BDI class. We found that this 
class also includes time-reversal invariant systems with a 
mirror symmetry (which was originally considered as a 
subset of the Dili class with mirror symmetry^). This 
includes for example all s-wave and d-wave superconduc¬ 
tors with spin confined to a plane and a homogeneous 
superconducting phase. In these systems it is possible to 
find an operator V which satisfies the anti-commutation 
relation {77, V} = 0, and to unambiguously identify a 
Majorana character, and thus two different types of Ma¬ 
jorana states. These are the eigenstates of V with eigen¬ 
values ±1. States with opposite character are localized 
and well separated. We have provided several examples 
of systems in which writing down the V operator is possi¬ 
ble, as well as the explicit form of the Majorana character 
operator for each example. It is also possible for the spin 
degree of freedom to be replaced by a sub-lattice degree 
of freedom, the same arguments then apply. 

For case (B) where no operator V satishes {77, V} = 0, 
it is impossible to ascribe a character to the Majorana 
states, these are the D and Dili topological supercon¬ 
ductors. We want to point out however that the labelling 
of the Majorana states, associated with their character, 
can be extended to systems where the properties hold 
at least locally, for example to SNS ring junctions along 
which the superconducting phase varies slowly, or to SNS 
junctions with different SC phases in the two SC regions. 
In as much as the conclusions of the character can still be 
applied there is a low energy effective theory in a higher 
symmetry class (BDI). 

We have also considered the consequences of the exis¬ 
tence of the two types on Majoranas on their hybridiza- 






tion/scattering and braiding properties. For case (A) it 
can be shown that states with opposite character can hy¬ 
bridize, destroying each other, whereas states with equiv¬ 
alent character can not hybridize provided {H,V} = 0 
holds. This is in agreement with the expectations from 
the topological invariant. Case (A) is always in the BDI 
class. In Id BDI has a Z invariant, allowing for mul¬ 
tiple Majorana bound states at a single end. For case 
(B) hybridization to a finite energy state is the norm, 
the exceptions being all finely tuned, though not impos¬ 
sible to find, examples. We note that case (B) is in fact 
the desired behaviour for braiding of Majorana states at 
junctions of wires. That the existence of multiple Majo¬ 
rana states is not possible is in agreement with the system 
being described by a Z 2 invariant. 

We have also considered the implications of case (B) on 
the braiding properties of the Majorana states by map¬ 
ping one characteristic {{H,V} yf 0) model to a Kitaev 
chain with inhomogeneous superconducting phase. We 
find that the non-abelian braiding remains intact. A 
good confirmation of this would be a numerical analy¬ 
sis of braiding on characterless Majorana bound states. 
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the ‘original’ basis of our problem, {ipi}- Then from 
{H, V} = 0 with V = cr^T^ we have V = VT resulting 
in 


V = (cos[0]cr^ — sin[0] cos[(^]cr’^ — sin[0] cos[(^]cr‘’) , 

(A2) 

acting locally in space (the spatial dependence will be 
suppressed where it is not explicitly required). The ro¬ 
tations in spin space, along with the change of phase be¬ 
tween the particle and hole spaces which follows, define 
the full 6-d space in which the set of possible V operators 
exist. 

An arbitrary but homogeneous superconducting phase 
of K requires the transformation V —>■ V G where 


/ e ‘«/2 0 0 0 \ 

0 e “/2 0 0 

0 0 0 

V 0 0 0 e-’''/^ y 


(A3) 


is the gauge transformation which sets the superconduct¬ 
ing phase to zero in the Hamiltonian. 

As we have a universal planar spin density then the 
angles 9 and (j) can be calculated from two, non-parallel, 
spin vectors at different spatial points. Let these be S'i_ 2 j 
then the spin density is confined to the plane perpendic¬ 
ular to the unit vector 


^1 X ^2 

1^1 X S2\ ’ 


(A4) 


Just prior to submission we became aware of similar re¬ 
sults on the nature of mirror symmetry in Ref. 38 derived 
by considering the Berry phase and the winding numbers 
of the topological phases, and to a review of topological 
classification and symmetries including spatial and mir¬ 
ror symmetries in Ref. 39. 

Appendix A: Planar spin configurations 

If we focus on a spinfull Bogoliubov-de Gennes (BdG) 
Hamiltonian with only real elements, (see for example the 
Hamiltonian in Eq. (7) where the superconducting phase 
is zero and the spin is confined to the ccz-plane), then 
V takes a particularly simple form: V = cr^r^. That V 
anti-commutes with the Hamiltonian follows immediately 
from the particle-hole symmetry and the realness of the 
Hamiltonian. In this case T+ — K and therefore V = 
C T+ is satisfied. 

By rotating from the trivial case we can find the generic 
operator V. We start with a system with some generic 
spin density which is confined to a plane. Gonsider a 
global rotation between two orthonormal bases ip = T'0 



We rotate the Majorana polarization from a basis in 
which the spin is confined to the ccz-plane, {ipi}, to 


and passing through the point (0,0,0). The necessary 
rotation is that which turns this into y, so that the spin 
density is confined to the xz plane and the Hamiltonian 
is real. Therefore we can write 

fis = (sin 6 cos (p, cos 9, sin 9 sin cp ), (A5) 

from which we can calculate 9 and (p. We can therefore 
rewrite the operator, Eq. (A2), as 

V = ^ , (A6) 

|5i X ^ 2 ! 

where naturally <t* refers to the complex conjugate of 
a. Global changes to the superconducting phase can be 
trivially incorporated, and we do not add them here ex¬ 
plicitly for the sake of clarity. 


Appendix B: Gauge transformation for magnetic 
inhomogeneity 

To understand why a spatially varying magnetic field 
can generate Majorana states it is helpful to make a 
gauge transformation and find the spin-orbit like terms 
which are generated by the gauge potential. 

This procedure is well known, but we repeat it here 






9 


for clarity and completeness. This gauge transforma¬ 
tion rotates the spin direction to that defined locally 
by the magnetic field. In such a basis the particles 
have an effective spin-orbit like coupling. A rotation 
^ja = 'Yha' ^ja' which diagonalizes the magnetic 

field term such that 

W] fij ■ ^Uj = (Bl) 


The new interesting term is generated in addition to Hq 
by the contributions like (^^U\ lAj Cj and is 


Hso = -l E (B5) 

<z,j'>,cr(7' 


is 



(B2) 


where 


This gauge transformation results in the Hamiltonian 

Ho + Hb + + H,^Ho + Hb+ Hso + Hs ■ (B3) 

The last term is generated by the gauge transformation 
from the hopping in Hq. The s-wave pairing term is in¬ 
variant under such a rotation. We note that d-wave pair¬ 
ing would be modified, a case we do not consider here. By 
construction Hb —t Hb, a diagonal homogeneous mag¬ 
netic field. More explicitly 

Ho + Hb+H, = Y, [(m - t)r^ + - Ar"] 

3 

H . (B4) 

<*j> 


S*-’ = icr*' sin cos - icr“ sin sin 

+\a-^ cos ipfj sin , (B6) 


and we have defined = {pi ± Pj)l‘2. and = 
{'di ± dj)/2. In 2-d systems a mapping between simple 
magnetic fields varying in the defined manner and spin- 
orbit coupling is not possible.^® For the 1-d systems we 
will focus on here it is straightforward. 
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